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Abstract 



The light cone QCD sum rules are derived for the K*Kix cou- 
q pling gx*K-K and the pmr coupling g p7T7T . The contribution from 

| 7 | ■ the excited states and the continuum is subtracted cleanly through 



the double Borel transform with respect to the two external mo- 
CnJ ! menta, pf, p\ — (p — q) 2 . Our result Qk'K-k = (8.7 ± 0.5) and 

g pnn = (11.5 ± 0.8) is in good agreement with the experimental 
value. 
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The pun (K*Ktt) coupling g p7T7T {gK*K-n) plays a very important role in the phenomeno- 
logical models for the nuclear force and nuclear matter. Now it is widely believed that 
QCD is the underlying theory of the strong interaction. Yet the complicated infrared be- 
havior of QCD causes the first principle derivation of hadron properties highly nontrivial. 
So a quantitative calculation of the g p7T7T (gK*Kn) coupling with a tractable and reliable 
theoretical approach proves valuable. 

The method of QCD sum rules (QSR), as proposed originally by Shifman, Vainshtein, 
and Zakharov [|TJ and adopted, or extended, by many others |2|, [| f|, are very useful in 
extracting the low-lying hadron masses and couplings. In the QCD sum rule approach the 
nonperturbative QCD effects are partly taken into account through various condensates 
in the nontrivial QCD vacuum. In this work we shall use the light cone QCD sum rules 
(LCQSR) to calculate the K*Kir and pirn couplings. 

The LCQSR is based on the OPE on the light cone, which is the expansion over 
the twists of the operators. The main contribution comes from the lowest twist opera- 
tor. Matrix elements of nonlocal operators sandwiched between a hadronic state and the 
vacuum defines the hadron wave functions. When the LCQSR is used to calculate the 
coupling constant, the double Borel transformation is always invoked so that the excited 
states and the continuum contribution can be treated quite nicely. Moreover, the final 
sum rule depends only on the value of the hadron wave function at a specific point, which 
is much better known than the whole wave function || . In the present case our sum rules 
invole with the pion wave function (PWF) (p n (uo = \)- Note this parameter is universal 
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in all processes at a given scale. In this respect, <p w (uo = \) is a fundamental quantity 
like the quark condensate. In || the value is estimated as <p n (uo = -) = 1.2 ± 0.3 using 
the pion nucleon coupling constant and the phenomenological pure coupling constant as 
inputs. Recently the light cone sum rule for g n NN{q 2 = 0) was reanalyzed JF). The 
contribution from the gluon condensate < g 2 s G 2 > and the quark gluon mixed condensate 
(g c qv • Gq) is added. The uncertainty due to A at is reduced in the numerical analysis with 
the help of the Ioffe's mass sum rule. A new value (^(1/2) = 1.5 ± 0.2 |7J is obtained 
using the experimentally precisely known g n NN @- 

The LCQSR has been widely used to derive the couplings of pions with heavy mesons 



in full QCD [fl, in the limit of mq — > oo [[| and 1/mg correction |T0[, the couplings 
of pions with heavy baryons |Tl[|, the p-decay widths of excited heavy mesons |12| and 
various semileptonic decays of heavy mesons f[3f etc. 

Our paper is organized as follows: Section [l] is an introduction. We introduce the 
two point function for the K*Ktt vertex in section |2|. The definitions of the PWFs are 
presented in section |3|. In the following section we present the LCQSR for the K*Kir and 
P7T7T couplings. A short summary is given in the last section. 

2 Two Point Correlation Function for the K*Ktt coupling 

The dominant decay mode is K* — > Kit for K* and p — > tttt for p. The relevant decay 
amplitudes are 

< K* (p)n-(q)\K-(p + q) >= -g K * K7T q^ (1) 

< K*-(p)n°(q)\K-(p + q) >= -^V* (2) 

< P'(p)^(q)K(P + q) >= -Qpnnq^ (3) 

where is the polarization vector of the vector meson. The resulting decay width reads 

r(K , ^ Kw) = (4) 

With T(K* -> Kir) = 50.8MeV and T(p -> tttt) = 154MeV |lj one gets g K , Kv = 9.08 
and g p7T7T = 12.16. 

To study these couplings, we start with the two-point correlation function: 

n(p) = i f (Pxe ipx {n-(q)\T[d(x)^s(x), s(0) 7 Q 7 5 w(0)]|0) . (6) 



The instanton contributions may invalidate the usual sum rule techniques for the pseudo- 
scalar current [|Tj|. We use the pseudo- vector currents for K and 7r throughout this work. 
At the phenomenological level the eq.(|6|) can be expressed as: 

with p\ = p, p2 = p + q. The ellipse denotes the continuum and the off-diagonal transition 
contribution. The decay constants fx and fx* are introduced as: 

(K( P )\s(0)W5u(0)\Q) = fKPu, (8) 
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(K*\d(0) ltl s(0)\0}=m K *f K ^ } (9) 

where e M is the polarization vector of the K* meson, 

After discarding the single pole terms in (0) which will always be eliminated through 
double Borel transformation later, the polarization operator can be expressed as 

nfoft, g) = g K *K. ( { K * mK ; u I KP \ M + ~(i - 4^1 • (io) 

(p 2 - m z K «) (jpt, - m z K ) 2 m z K * 

In the following we shall focus on the tensor structure p a q ti . 

3 The Formalism of LCQSR and Pion Wave Functions 

Neglecting the four particle component of the pion wave function, the expression for 
H(px,p2,q) with the tensor structure at the quark level reads, 

U(p 1 ,p 2 ,q) = -iJ rf 4 xe^Tr{(7r-(g)| M (0)J(x)|0) 7 ^^(x)7 Q 75} • (11) 

where iS s (x) is the full strange quark propagator with both perturbative term and con- 
tribution from vacuum fields. 

d 4 k k + m s 



iS s (x) = (0|T[s(x),s(0)]|0) = -i / -—e 



(2tt) 4 (m 2 s - k 2 ) 
( ' S ' S) -^(s9svGs) + --- . (12) 
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where we have introduced k = fc M 7 M , m s = 150MeV is the strange quark mass, _D M = 
- ig s A^. 

By the operator expansion on the light-cone the matrix element of the nonlocal op- 
erators between the vacuum and pion state defines the two particle pion wave function. 
Up to twist four the Dirac components of this wave function can be written as jSf: 

< 7^)1^)7^(0)10 > = -iUq, f 1 du e™«%» + x 2 9l (u) + 0(x 4 )) 

J 

+ U(x, - ^) C du e iu « x g 2 (u) , (13) 
qx Jo 

— f to 2 

< n(q)\d(x)ij 5 u(0)\0 > = Jw n / du e iuqx y P (u) , (14) 



m u + m d Jo 

— f m 2 

<-K(q)\d(x)a^ 5 u(0)\0> = %(q^x v -q v x ll ) . U * du e iwpi <p (u) . (15) 

6(m u + m d ) Jo 



< ■K(q)\d(x)a a pj 5 g s G tlv (ux)u(0)\0 >= 

iM(qi*Qa9vi3 - q v qa9rf) - (q^gua - q^qpg^)] J ^(ai)e iqx{ai+va3) , (16) 

< 7r(g)|J(x)7 A1 755f s G Q/3 (fx)M(0)|0 >= 
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, r / ^ a qfj,\ t 

U [9/3 [pop - ^7 ) - ?a , 

q ■ x J 



(17) 



and 



< 7r(g)|d(x)7^G a/3 (ux)M(0)|0 >= 



/vr[g/3( 



) -q a [9fr ) 

q ■ x ' K q ■ x ' 



q ■ x 



(q a xp - q p x a ) / Va^a^e^^ . 



;is) 



The operator G a p is the dual of G a p\ G a/ 3 = ^ a /3SpG 5p ; T^ai is defined as T>oti = 
da\da 2 da^5{l — a\ — a 2 — 03). Due to the choice of the gauge x^A^x) = 0, the path- 
ordered gauge factor Pexp (ig s dux^A^ux)) has been omitted. The coefficient in front 
of the r.h.s. of eqs. (0), ( p~5|) can be written in terms of the light quark condensate 



<uu > using the PCAC relation: fi„ 



ml 



-y i <UU>. 

■1 7T 



m u + m d 

The PWF (p n (u) is associated with the leading twist two operator, gi(u) and g 2 {u) 
correspond to twist four operators, and <pp(u) and (p a (u) to twist three ones. The function 
V^s-jr is of twist three, while all the PWFs appearing in eqs.fll7|), fll8D are of twist four. The 
PWFs (p(xi,fi) (n is the renormalization point) describe the distribution in longitudinal 
momenta inside the pion, the parameters x« — 1) representing the fractions of the 

longitudinal momentum carried by the quark, the antiquark and gluon. 

The wave function normalizations immediately follow from the definitions (|13|)-(|T8|): 
fodu ip n {u) = fodu ip a (u) = 1, Jodu gi{u) = S 2 /12, JVa^^ai) = J Vaiip\\{aii) = 
0, / Va^xXai) = — J Uati<p\\(ai) = P/3, with the parameter 6 defined by the matrix 
element: < n(q)\dg s G a ^ a u\0 >= i^j^. 



4 The LCQSR for the K*Kn coupling 

Expressing ([□]) with the strange quark propagator and keeping only the tensor structure 
p a q fl , we arrive at: 

1-1 du 



f du 

n(pi,P2,g) = -i/w / — n — 7 — ; — ^{vA u 

Jo mi — (v + uar 



?f^ a {u) + Aug 2 {u) - A gi {u) - AG 2 {u) 8m 2 s [ gi (u) + G 2 (u)} 



m 2 s — (p + uq)' 

+ ■■■} 



m 2 — (p + uq) 2 \m 2 s — (p + uq) 2 } 2 

where /j, n = 1.65GeV, f m = 132MeV, G 2 {u) = — J M g 2 (u)du, which arises from integration 



by parts to absorb the factor l/(q ■ x), 

-1 „— iuq-x r \ 



q ■ x 



g2 (u)du = -i I e- iuq - x G 2 (u)du - G 2 {u)e- iuq x \l , (19) 



Note the second term in (|19|) vanishes due to G 2 (uq) = at end points uq = 0, 1. 

Making double Borel transformation with the variables p\ and p\ the single-pole terms 
in (|7]) are eliminated. The formula reads: 
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Subtracting the continuum contribution which is modeled by the dispersion integral 
in the region s 1? s 2 > s 0; we arrive at: 

2,2 22 

m K" t " m K* '"q( 1 ~ u o)y + m s 9 _ _fl 

f K fK*m K *g K * K7T e ^ = f n e {M (1 - e ^^(tto) 

+ ^/i 7r <^M + 4u ^K) - 4^i(w ) - 4G 2 (n ) - -^#[3iK) + G 2 (u )]} , (21) 

where s\ is the continuum threshold, u = M ^ M -i , M 2 = ^^m 7 ^ Mf, M% are the Borel 
parameters. Note all the PWFs involved with the vacuum gluon field do not contribute 
to the tensor structure p a q 13 , which greatly simplifies the analysis of the sum rule. 

Similarly we can obtain the LCQSR for pirn coupling. In our calculation we take the 
up and down quark current mass to be zero. Then we arrive at the simple sum rule: 

f P m p g pn7T e w = V2e m 1 {M (1 - e i 7 )^^) 

+Au g 2 {u ) - A gi {u ) - AG 2 {u )} . (22) 



The sum rule (pip and (|22|) appears asymmetric with the different masses for K and 
X* mesons at first sight. However, if a sum rule holds well, there should exist a certain 
interval called the working region of the Borel parameter, Mf < M\ < Mj. Within this 
region the sum rule is insensitive to M\ and stays reasonably stable with the variation 
of Mb- In other words, every point of Mb in the above interval is equally good for the 
analysis of the sum rule. So long as the working regions for M\ and Mf have some 
overlapping region, which does occur in our case, we can choose a common value in the 
overlapping region for both Mf and Mf . From the above argument we know such a 
choice will not alter significantly the final result of the sum rule, i.e., the choice Mf = Mf 
is allowed in the analysis of the LCQSR. Moreover, the symmetric choice of Mf = Mf 
will enable the clean subtraction of the continuum contribution, which is crucial for the 
numerical analysis of the sum rules. In contrast, the asymmetric choice will lead to the 
continuum subtraction extremely difficult in ||, which is the operational and technical 
motivation for the choice of Uq — 1/2. We shall work in the physical limit q 2 = m 2 — > 
in (0). 

The resulting sum rule depends on the PWFs and the integrals of them at the point 
u = |. We adopt ip n (u ) = 1.5 ± 0.2 0. For the other PWFs we use the results given in 
§, ip^uo) = 1.47, gi(uo) = 0.022GeV 2 , g 2 {u ) = and G 2 {u ) = 0.02GeV 2 at u = \ at 
the scale fi = lGeV. 

The overlap amplitudes fx* and f p can be determined in a self-consistent manner 
making use of the corresponding mass sum rules. For example, for the p meson 0,0], we 
have: 

^ 2 ^% + >^-s + v^i>' < 23 > 

where E\ = 1 — (1 + j$r)e m b (with the continuum threshold s = 1.5GeV 2 ) is the 

factor used to subtract the continuum contribution. Numerically, f p = (0.18 ± 0.02)GeV, 
f K . = (0.21 ± 0.02)GeV, f K = (0.15 ± 0.02)GeV |, @, 0, |T§/ 

Note that we have used the pseudo- vector interpolating current, which couples strongly 
to both pseudo-scalar mesons n, K and pseudo- vector mesons ai(1260), i^i(1270). ai(1260) 
is a broad resonance with a full width of ~ 200MeV. In order to eliminate the contami- 
nation from ai(1260) and fTi(1270), we choose the continuum threshold parameter to be 
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s 2 < (m ai - ^) 2 ~ (1.2 ± 0.1)GeV 2 and sj < (m Kl - ^) 2 ~ (1.5 ± 0.1)GeV 2 These 
values are consistent with the continuum threshold for the sum rules of pseudo-scalar 
mesons. But they are slightly smaller than the continuum threshold for vector mesons. 

The final sum rules are stable with reasonable variation of the Borel mass around m? K * 
and m 2 respectively after the exponential factor is moved to the left hand side as can be 
seen from Fig. 1 and Fig. 2. The terms with the strange quark mass m s contribute 
about 10% to the whole sum rule (pip. Our final result is Qk'K-k = (8.7 ± 0.5) and 
g pn7T = (11.5 ± 0.8) with the central values of f p , fx, fx*, which agrees very well with the 
value extracted from the experimental data gK*Kix = 9.08 and g p7T7T = 12.16 |14 . 



Although the coupling of K* pK may be sizeable, the pK intermediate states con- 
tribute to the continuum only because m p + m K = 1.27GeV > y/si. In our approach we 
have invoked quark-hadron duality and modeled the spectral density with s > si at the 
hadronic side with the free parton-like one. So the pK contribution is subtracted away. 
On the other hand, the intermediate state 2>nK does lie around itik* and below the con- 
tinuum threshold y/s{. But due to the strong suppression from the four-body phase space 
integral, its contribution is negligible. In other words, the K* pole term dominates the 
possible intermediate states with same quantum numbers below the continuum threshold 
Si while those intermediate states above S\ is subtracted away using the quark-hadron 
duality assumption, which is the corner stone of the QCD sum rules approach. 

In summary we have calculated K*Kir and pmr couplings using the light cone QCD 
sum rules. These couplings are related to the values of PWFs at the point uq = \, which 
are universal in all processes. Our results are in good agreement with the experimental 
data. 
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Figure Captions 

FIG 1. The sum rule for g p ^^ as a functions of the Borel parameter M 2 . From bottom to 
top the curves correspond to S2 = 1.1, 1.2, 1.3GeV 2 . 

FIG 2. The sum rule for qk*k-k as a functions of the Borel parameter M 2 . From bottom 
to top the curves correspond to si = 1.4, 1.5, 1.6GeV 2 . 
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